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^^ ' We show that the statistics of tunnelling can be dramatically affected by scarring and 

Cn , derive distributions quantifying this effect. Strong deviations from the prediction of random 

matrix theory can be explained quantitatively by modifying the Gaussian distribution which 

describes wavefunction statistics. The modified distribution depends on classical parameters 

which are determined completely by linearised dynamics around a periodic orbit. This 

distribution generalises the scarring theory of Kaplan [Phys. Rev. Lett. 80, 2582 (f998)] 

to describe the statistics of the components of the wavefunction in a complete basis, rather 

than overlaps with single Gaussian wavepackets. In particular it is shown that correlations 

in the components of the wavefunction are present, which can strongly influence tunnelling- 

\^ \ rate statistics. The resulting distribution for tunnelling rates is tested successfully on a 

(-H ■ two-dimensional double-well potential. 

I. INTRODUCTION 

^- ' Tunnelling rates from systems with complex internal behaviour are naturally described statistically. Such 

analyses have been used in particular to understand nuclear resonances [|J , chemical reaction rates [g-|| and 
quantum dots P,E0[. While the manner in which internal states couple to the continuum varies in these 
analyses, the internal quantum mechanics is usually modelled using one of the standard ensembles of random 
matrix theory (RMT). We show here that in systems where the complexity of the internal dynamics derives 
from low-dimensional chaos, sufficiently strong deviation from the standard RMT statistics is possible that 
it dominates the statistics of tunnelling. These deviations were pointed out in pT] a nd are explained in 
detail here using the theory of scarring developed by Heller, Kaplan and coworkers |12|]l3| . 

We consider tunnelling to and from regions of phase space associated with chaotic behaviour. In particular, 
our theory works for the calculation of level splittings in double well potentials and of resonance widths 
of quasibound states in metastable wells. With the development of experimental techniques capable of 
measuring state-specific reaction rates, such tunnelling from highly-excited chaotic states has become directly 
relevant to the analysis of chemical reactions, for example. Specific examples of reactions which have been 



> 



x 

treated using RMT-based statistical methods include the dissociation rates of D2CO |2|,|3| and N2O ]8|]. 
In |11| , the distributions arising from such an analysis were shown to be determined in a simple way by 
the stability and action of a complex tunnelling orbit which crosses the potential barrier with minimum 
imaginary action. The resulting statistical distributions for the tunnelling rate agree well with numerically 
computed ensembles except when the real extension of the tunnelling orbit into the potential well is periodic; 
in that case, strong deviations from the RMT prediction are observed and it was proposed in Q that these 
are due to the effect of scarring on wavefunction statistics as outlined in |13| | . Additional evidence in support 
of this has subsequently been provided in |1J] . 

In this paper we derive a distribution describing scar-influenced tunnelling-rate statistics in two- 
dimensional potentials. It is determined completely by three dynamical parameters which are calculated 
from the monodromy matrices of the tunnelling orbit and of its real periodic extension. These three param- 
eters can be understood as the stability of the scarred orbit, the stability of the complex tunnelling orbit 
and an angle relating stable and unstable manifolds. The root of this calculation is a conjecture governing 
the statistics of the components of chaotic wavefunctions in the eigenbasis of a tunnelling operator which 
was defined in |laj . The low- lying states of this tunnelling operator are approximated by the eigenstates of 



a harmonic oscillator. Following the linear scarring theory of Kaplan |13| we show that the components of 
chaotic states in this basis have variances which deviate from those of RMT and which must furthermore be 
correlated. We conjecture that the eigenfunction distribution is a nonisotropic Gaussian. This is the simplest 
distribution consistent with the observed correlations and may also be obtained from maximum-entropy ar- 
guments. This conjecture is tested against quantum-map models and found to describe statistics of their 
eigenstates very accurately. The conjecture is then used to derive modified distributions for tunnelling rates 
in the presence of scarring and these are found to describe well the statistics of energy-level splittings in 
chaotic double well potentials. 

We see the benefits of this work as being twofold. First, there is a solution to the primary problem of 
incorporating the effects of scarring into the statistics of tunnelling. This enables the detection and interpre- 
tation of system-specific dynamical detail in measurements of multidimensional tunnelling and furthermore 
provides a very direct manifestation of the phenomenon of scarring in quantities with real physical relevance. 
Second, we believe that from the point of view of scarring alone, even without reference to tunnelling, the 
correlated joint probability distribution proposed here for the statistics of eigenstates holds considerable 
promise as a theoretical tool. In particular, it suggests methods of quantifying the effects of scarring in a 
way that does not depend on a choice of test states. This second aspect is not fully exploited in the present 
paper, but we believe that the essential elements needed for a full development are put in place. We remark 
that the effect of correlations on the statistics of tunnelling rates was also considered in a somewhat different 
context in p0[ . 

An outline of the paper is as follows. We begin in section || with a review of the existing theory of the 
statistics of tunnelling rates. This includes a definition of the tunnelling operator introduced in [ p.5[ and a 
summary of the derivation from it of nonscarred tunnelling-rate distributions as outlined previously in |Ll| . 



In section III we use the linear theory of scarring to deduce that in the presence of scarring the components 
of chaotic eigenstates in an eigenbasis of the tunnelling operator must be nonRMT and conjecture a modified 
Gaussian distribution for them. This conjecture is tested on quantisations of perturbed cat maps and found 



to work well. In section IV we use the conjecture to deduce modified tunnelling-rate distributions in the 
presence of scarring and these are compared successfully with distributions of splitting in chaotic double- well 
potentials. 

II. TUNNELLING STATISTICS WITHOUT SCARRING 

We begin with a brief description of existing theory of tunnelling-rate statistics in the absence of scarring, 
on which our present work is based. We first describe a tunnelling operator defined in |jj| which relates 
tunnelling rates to the properties of wavefunctions in the allowed region. We then outline how this tunnelling 
operator is used to derive distributions for tunnelling rates in the absence of scarring. 

A. The tunnelling operator 

A full description of the semiclassical matrix element and complete definitions of the various objects needed 
to calculate it are described in fl5|| . For the purposes of statistical analysis, an idealised model formulated 
in terms of quantum maps suffices, which we summarise here. 

Our model starts with a finite-dimensional Hilbert space Ti. We take this space to be the quantum analog of 
a Poincare section E in a potential well []16[ , with the dimension N of H being approximately proportional 
to the area of E. A quantum mapping acts on the space Ti as a unitary operator U, whose eigensolutions 
we denote by 

U\n) =e- w "\n) 



and whose classical analog is a symplectic map F corresponding to chaotic motion within the well. In |15|| , 
the eigenstates \n) are represented as certain cross-sections of the eigenfunctions of the usual Hamiltonian 
in the full Hilbert space. 

Tunnelling rates are calculated using a tunnelling operator T which also acts on Ti and from which the 
scaled tunnelling rate associated with a chaotic eigenstate \n) is calculated from 

Vn = ~^ r - (1) 

In this formula, we normalise the chaotic eigenstates according to (n\n) — N so that, in particular, the 
overlaps with a second basis |fc) of Ti satisfy (|(fc|n)| 2 ) = 1. As a result, (y n ) = 1. In resonance problems, 
this scaled tunnelling rate is y n = T„/ '{T n ) where T n is the resonance width associated with a particular 
metastable state labeled by n and in double- well problems y n = AE n / (AE n ) where AE n is the splitting of 
a doublet labeled by n. 

The tunnelling operator T is interpreted as an evolution operator whose classical correspondent is a 
symplectic map T which is complex. This map T has a real fixed point Co which corresponds to the most 
probable tunnelling route. That is, using Co as an initial condition, a complex trajectory of imaginary period 
can be traced out which crosses a potential barrier with minimum imaginary action. Mappings of points 
in a neighbourhood of Co i n the Poincare section E correspond to complex trajectories in full phase space 
which cross the potential barrier near this central tunnelling route. The matrix element in (fy) then samples 
the state \n) in a small region of S surrounding Co and with area of 0(h) (as determined by say, a Wigner 
function on E |lq]). This region is determined by the dynamics of T near Co and outside of it the symbol 
of T decays exponentially. For this reason it is consistent within semiclassical approximation to linearise 
dynamics around Co and let T be the quantisation of the resulting complex symplectic matrix W. Note that 
in practical terms W is simply the monodromy matrix of the tunnelling orbit. 

We now describe the tunnelling operator for resonance problems of symmetric double wells where the 
symmetry is a reflection (x,y) i—*- (—x,y). If we choose a representation in which Co is at the centre of 
coordinates C = {q,p) on E, we may write 

q- _ —aoh/h 

for these cases, where h is quadratic in the corresponding operators C = (<?>?>) and a > . The corresponding 
symbol h(q,p) is a quadratic form in {q,p) 1 expressed as 

H°,p) = \(" K ^ 

where K is a positive-definite symmetric matrix. The matrix K is determined from W by writing 

W = e ~ ia ° JK . 

where J is the unit symplectic matrix. We normalise K so that det K = 1 and then the parameter ao is 
fixed by a calculation of W. That W can be written in this way for real ao and K can be shown using the 
conjugate-time-reversal symmetry 

W* = W~\ 

which is a special case of a similar symmetry of T [ Jl5[ . This symmetry is a reflection of the fact that the 
tunnelling orbit is periodic with an imaginary period and so complex conjugation is equivalent to time- 
reversal. Note that this is the case even if the problem does not have a proper time-reversal symmetry (such 
as in the presence of magnetic fields, for example). 

We limit the detailed discussion in this paper to two-dimensional potentials, so that the Poincare section 
E has one degree of freedom. In that case we may denote the eigensolutions of h by 



h\k)= (k + ~\h\k), 

using a tilde to distinguish them from the chaotic eigenstates \n). Note then that the eigensolutions of T 
are then 

where A = e a ° . Notice that A is the larger of the two eigenvalues e ±Q0 of W (which are real even though 
W is complex). Finally, we note that in double well problems for which the symmetry is an inversion 
(x, y) i— ► (—x, —y), the eigenvalues of W as defined in jl5[ are negative and the eigensolutions of the tunnelling 
operator are of a similar form except that the eigenvalues are A~ fc |A| ' . 

B. Statistics using standard RMT 

A standard statistical model of eigenstates of chaotic systems is that their components in a generic basis are 
Gaussian-distributed. This assumption has been combined with various models of coupling to the continuum 
[p 10 O] to produce statistical distributions for the scaled tunnelling rates in the absence of scarring. The 



discussion in [ p.l| forms the basis for our treatment of scarred distributions and we therefore summarise it 
here. We assume for ease of presentation that the system is time-reversal symmetric so the assumptions of 
the GOE are adopted. Results for GUE systems are derived similarly and are summarised at the end. 
We expand the eigenstate \n) in an eigenbasis of T, 



l ) = Yl x k\k), 



suppressing the dependence of the coefficients x& on the chaotic-state index n. The normalisation of \n) is 
such that 



|Xfc| 



and in systems with time-reversal symmetry, the Xfc's are real. The standard RMT model is that for large 
N, the Xfc's are uncorrelated and distributed with a joint probability distribution 
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k=0 



J_ e-^/ 2 
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1 e- xTx / 2 , (2) 



where x = (xo, xi, X2, • • •). Given that (gj) expresses each y as a quadratic form 

y = x T Tx 

in x, where T is an N X N matrix representing T/TrT in a |fc)-basis, we can write the distribution for y in 
the form 

P(V) = J S(y- x T Tx) P(x) dx. (3) 

Fourier-transforming with respect to y gives the characteristic function, 



p(q)= / e i "yp(y)dy 



3C: 
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exp 



(2n) N / 2 
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x T (I-2iqT)x 
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dx 



A/det (/ - 2i g T) 
One finds similarly that 



(4) 



^ = det (J -igT) (5) 

in the absence of time-reversal symmetry. 

These distributions have been shown in 1 1 to describe successfully the tunnelling rate statistics of chaotic 
double wells when the tunnelling route has a nonperiodic real extension. When the tunnelling route connects 
to a periodic orbit, however, strong deviations are found [ |ll| , [l4| which we explain in the next sections using 
the idea of scarring. 

III. WAVEFUNCTION STATISTICS AND SCARRING 

It has been pointed out by Kaplan and coworkers that the Gaussian distribution in (g) does not describe 
wavefunction statistics if a basis state \k) is localised in phase space near a periodic orbit of the real classical 
dynamics (corresponding to the Poincare map F). In tunnclling-rate statistics, this can lead to strong 
deviations from the distribution described the previous section if Coi in addition to being a fixed point 
of the complex map J-", is a fixed point or short periodic orbit of F. In this section we review elements 
of that argument and use them to calculate correlations among the components of x. From this analysis 
there emerges a natural conjecture for generalising (g). The distributions we calculate on the basis of this 
conjecture are necessary to understand tunnelling but also, independently of the tunnelling problem, hold 
promise as a basis for understanding the statistics of scarring in a rather general way. 

A. A Gaussian hypothesis 

For a given ensemble of states with overlaps x = (xq, X\, X%, ■ ■ •), let us define a matrix of correlations C 
whose elements are 

Cik = (xix* k ). (6) 

The average is over the eigenstates \n) and will be made more precise below. We will in particular consider 
ensembles formed by the eigenstates \n) of U whose eigenangles lie in a subset of the unit circle. By varying 
the quantum dimension N, for example, we can consider ensembles corresponding to fixed classical maps F 
and T . If the joint probability distribution is of the form (g) or its GUE equivalent, we expect 

C = I, (7) 

which says simply that (|xfc| 2 ) = 1 and that Xk and xi are uncorrelated unless k = I. In order to avoid 
confusion with a correlation function defined below, we will henceforth refer to C as the covariance matrix. 
We now show that deviations of C from the identity can be calculated on the basis of linearised dynamics 
around the point Co- The construction begins with a calculation of quantum recurrence for basis states 
localised near Co and is finished by a Fourier transformation which relates these to the averages {xix* k ). 



Let {|fc)| k = 0, 1, • • •} be an eigenbasis of T as defined in the previous section and suppose that Co is a 
fixed point of the classical map F corresponding to U. We now argue that for some ensembles (R) cannot 
hold. Following |1J] we define the correlation function 



A lk (t) = (l\U f \k) 

N 



^5>aSe-""« (8) 

n=l 



and denote its Fourier transform by 



Si k (0)= J2 ei6tA ik(t) 



<2ir N 

= jr52*i*iW-9 n ). (9) 

n=l 

The correlation function is approximated semiclassically by a sum of the form 



over classical orbits p, typically complex, which begin in and return to a small region around Co after t 
iterations of F, assuming k and I are 0(1). Each orbit contributes a term with an amplitude A p which 
is independent of N and an exponential exp[2w NiS p ] with a rapidly-varying complex phase which decays 
exponentially with the distance of the contributing orbit from Co- The primary such orbit corresponds to 
the fixed point Co itself. We are free to choose the phase of U so that the action of Co vanishes (this amounts 
to a choice of the origin of the 9 n values). With that choice it is manifest that its contribution to Au-{t) is 
iV-independcnt . 

Consider averaging Aik(t) over a sequence of quantum systems with varying N but identical classical 
limits. We find then that 



(4*(t)>« y> P (e 2 ^ Sp ) 



and, with the exception of the fixed point Co f° r which S p — 0, the contributions average to zero. The 
contribution from Co can be obtained from linearised dynamics and we find that 

(A lk {t)) a A^(t) = (l\UU~k) 

where Un n is the quantisation of a linearisation of the real map F at Co- The linearised correlation function 
Af^(t) is calculated in the next section. It is iV-independent and decays exponentially to zero with time 
when k and I are 0(1). 

Performing the same averaging procedure on Sik (0) gives 

oo 

(S lk (9)) « S\m = E e-^A^W. 

t — — oc 

This is a real, 27r-periodic function of 9 which is also iV-independent. Suppose we form an ensemble of 
chaotic states by varying N and selecting states for which 9 n lies within a narrow window centred about 9. 
The eigenfunction expansion in (J9j) indicates that 



C lk {9) = ( Xl x%) = (S lk (9)) » S£(0). (11) 

Wc find therefore that, if we produce an ensemble of overlaps from states within a narrow window of phase 
angle, the covariance matrix violates condition (M) and the joint probability distribution for x necessarily 
deviates from the standard RMT form. 

The simplest distribution consistent with (Thy) is the Gaussian 

P(x;0) = m e -/3x-c-VW2 (12) 

where /3 = 1 in the GOE case and (3 = 2 in the case of GUE. We conjecture that these distributions 
govern the wavefunction statistics in the presence of scarring when states are taken from a small window 
of eigenphases. A full justification of this would require a detailed analysis of long orbits contributing to 
( |i"0| ) , along the lines of the nonlinear theory of scarring in [O . We will not offer such an analysis here and 
assume dl3) on the grounds of simplicity alone. This may be formalised by noting that (J13) minimises the 
information content [|l7|,|l8l 



X[P(x)] = /p(x)lnP(x) dx, 



subject to the constraint that the correlations are as written in (|g) and that N is large. Minimisation of 
information is not a proof, however, and ultimately we rely on detailed numerical testing to justify the choice 
of distribution. We now show in detail how C{9) may be calculated and verify that the Gaussian hypothesis 
describes quite well the scarred wavefunction statistics of specific quantum-map models. 

B. Calculation of the covariance matrix 

Construction of the scarred distribution in ( |12|) begins with the calculation of the linearised correlation 
function Afj^(t), from which the covariance matrix C{9) is obtained by Fourier transformation. We show in 
this subsection how A|™(£) may be calculated from the monodromy matrices linearising real and complex 
dynamics around the fixed point £ (corresponding to the maps F and T respectively) . In the main text we 
present enough detail that the method of calculation of Af^(t) should be clear and leave the derivation and 
other detailed discussion to appendices. Once A\™(t) has been constructed as outlined below, the covariance 
matrix is then easily calculated in practice using a FFT, for example. 

It is shown in appendix M that A\^(t) vanishes if I and k are not both even or both odd and may otherwise 
be calculated from the polar form 

4J+2„,fcW = G kn (ljj(t)) C <(*+V2)(*(t)-Mt*)44n*(t)+W2 j (13) 

where the Maslov index /i and the angles ip{t) and <j>(t) are defined by dynamics around Co as described 
below. The amplitude is 



where Cj? (x) denotes a Gegenbauer polynomial. In the special case n = we get the autocorrelation function 



All(t) = ^/cosi{;(t) P k (cosi;(t)) e < (*+ 1 /a)(*(t)-Mfar) j 



where P k {x) is a Legendre polynomial. We assume n > and t > in ( [13|) and use Af k n (— t) = Af£(t) = 
A l ^(t) to calculate the correlation function when t < or / < k. Note that this assumes a convention for 



the phases of the eigenstates |fc) that is outlined at the end of appendix H. In particular, this convention 
leads to a covariance matrix C(0) which is real-symmetric even in the case of GUE statistics; that is, while 
individual values of Xix* k are complex, the phase convention is such that the averages are real. 

To define the angles (/>(£) and ip(t), let Mo be the real symplectic matrix linearising the unstable real 
dynamics around Co and let 

W{z) = e l{lnz)JK (14) 

be a generalisation of the complex symplectic matrix described in the the last section, reducing to that case 
when z = e~ Q °. Euler expansion of the exponential in (JT3) leads to the identity 



TrW{z)M t = (-1)^ (m(t)z + 



m*(t) 



where 

m(t) = cosh pt + iQsmhpt 

and where (— l) A 'e ±p are the eigenvalues of Mq. The real parameter Q follows the notation of |13| and is 
given by Q = cot ip, where ip is the angle between the stable and unstable manifolds of Mo- Note that the 
parameter Q depends implicitly on W because the angle ip is calculated using the metric defined on phase 



space by K. Finally, we complete the explanation of (13) by defining the angles tj){t) and <p(t) implicitly 
using the polar decomposition 

m(t) = secip(t)e i<Kt) 

oim(t). We may choose these angles to lie in the ranges < ip(t) < ir/2 and —i^it) < <f>(t) < i/j(t) respectively 
(note that sec iji cos <f> — cosh pi > 1) and, with this restriction on <j>{t), the half-angle in (|13| ) is well-defined. 
Note that, as t — > oo, we find 

2 
cosib(t) e~ pt 

and the correlation function in ( |l3|) therefore decays at the exponential rate e~ pt ^ 2 when k is even and at 
the rate e _3p */ 2 when k is odd (since C k is then an odd polynomial). Less obvious from (|l|) is the fact 

that the correlation function also decays for fixed t with increasing k and I — k + 2n. This is a reflection of 
the fact that more excited states of the oscillator h are less localised at the fixed point Co and therefore less 
affected by scarring. 

When the linear correlation function is Fourier-transformed, we therefore get a covariance matrix Cik{0) 
which can be decomposed into two blocks corresponding to odd and even I and k and whose elements 
approach those of the identity matrix (Cik(9) ~ + Sik) as I and k increase. We expect scarring therefore 
to affect the statistics of Xk significantly for relatively small values of k and to have little effect when k is 
large. This exponential convergence towards identity means that the matrix C — I is trace-class and the 
determinant of C exists in the limit N — > oo. In the next subsection we calculate the resulting covariance 
matrix for some model quantum maps and show that the Gaussian hypothesis provides a good model for 
the wavefunction statistics of scarred systems. 

C. Testing the Gaussian hypothesis 

We test the Gaussian hypothesis using a model system in which U quantises the perturbed cat map, 



F e = F oM e , 



where 



is the standard cat map and 



1 1 
1 2 



mod 1 



M 



mod 1 



q \ / q — e sin 27rj? \ 

.p/ V p / 

is a kicked Harper map. For sufficiently small values of the parameter e this map shares the hyperbolic 
structure of the unperturbed map but not its nongeneric degeneracies. Details of the quantisation of these 
maps for Hilbert spaces of arbitrary dimension N may be found in |19| , [20| . 

The map F e has a symmetry of inversion about the origin, which is always therefore a fixed point. It also 
has a less obvious time-reversal symmetry. One finds that 



PF f p = f: 



where 



P 



-1 
-1 



mod 1. 



Note that P is antisymplectic and that P 2 = I. The map therefore allows us to test statistics in the GOE 
scenario. (A systematic discussion of such time- reversing symmetries can be found in pl[ ). 

To complete our model of a tunnelling system, we construct an analog of the tunnelling operator. Consider 
the Harper-like Hamiltonian 

h(q, p) = 2 — cos 7r(2p — q) — cos 27rg. 

This is invariant under the time-reversal symmetry P and has a minimum at the origin in phase space. The 
function does not have the full periodicity of the torus in q but, by placing the resulting discontinuity at 
q = 1/2, may be quantised in such a way that the ground and first few excited states of its quantisation h are 
localised near the origin. These are then approximated semiclassically by those of the harmonic oscillator 
(ft ~ <?/2) 2 + q 2 ■ The statistics of the components Xk of chaotic eigenstates of U in an eigenbasis of h therefore 
provide a test of the Gaussian hypothesis outlined in the previous subsection and since U and h share a 
time-reversal symmetry, the GOE hypothesis is appropriate. 
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FIG. 1. The smooth curves are calculated using our prediction for the covariance matrix as a function of phase 
angle 9. We show Coo(d), Cq2(6) and C'22(9) in (a), (6) and (c) respectively. In each case we find good agreement 
with the step curve, which represents the local average obtained by binning components of the eigenstates of U for 
e — 0.1 and N between 100 and 700. The points represent overlaps of individual states. 



We first consider the prediction for the covariance matrix C(9). The element C 00(d), which corresponds 
to the overlap statistics of a Gaussian wavepacket (defined as the ground state of h), has been considered 
extensively by Kaplan and coworkers. For comparison with later calculation it is plotted as a function of 9 
in Fig. 0(a) for the map F e with e = 0.1. Following Jl3[ , we refer to this curve as the spectral envelope. We 
expect larger than average overlaps where Coo(0) > 1 and we will refer to this as the scarred region. The 
rest of the interval is the antiscarred region, where we expect smaller than average overlaps. The points in 
the figure give a representative sample of individual overlaps. There is also a curve representing a binned 
average of these overlaps which is obtained by accumulating quantum data for even values of N between 
100 and 700. It is also important for our analysis to understand the statistics of xix k when / and k lable 
excited states of h. We therefore make an analogous comparison with quantum data for Co2(9) and Cfi(9) 
in Figs. Wb) and p](c) respectively. Note the different vertical scales in these cases. In all cases the average 
obtained by binning quantum data agrees well with the prediction Cik(0) obtained by Fourier-transforming 
the correlation function. This is to be expected since we have a relatively direct argument showing that this 
should be so. Less obvious is that the distribution of individual overlaps about these averages should be 
governed by the Gaussian hypothesis, since we do not have a proof for that. That this assumption can be 
verified numerically (below) is therefore nontrivial. 

To compare the distributions about the mean Ci k (9), lets us denote 



£ik = Rexix* k = xT5; fe x, 

where Bi k is a symmetric matrix whose ij th entry is (5u5j k + Si k 5ji)/2. (For the moment we allow either GUE 
or GOE statistics even though all the numerical calculations we present are for the GOE case). Similarly to 
the calculation of p(y) in section O, we find that the local distribution p(£i k ; 9) is predicted by the Gaussian 
hypothesis to have the Fourier transform 



Pik(q;d)= e^'"p^ lk ;9)d^ k 

J — OO 

1 

" [(2ir/f3) N detC(0)f /2 



.lot ( / - ^-B lk C(9) 



exp 

-0/2 



-\*\c-\6)- 2 fB lk ), 



dx 



(15) 



The determinant simplifies on expansion to 



2iq 



det ( I - -j-B lk C(9)\ = 1 - -±Ci k 



2iq 



-Tp (CuC kk -Cf k ) 



and Fourier transformation of (|15|) then gives 



p(Zu.;9) = e^" c ^ D ^f(^ lk ;D lk ,B lk ). 



Here, D, 



Ik 



CuC kk — C lk , B\ k 



\JCuC kk and /(£; D, B) is an even function of £, defined by 

1 

m-,D,B) 



nD 1 / 2 
1 



K (B£/D), for GOE and 



B 

(in which Kq(z) is a modified Bessel function). 



exp [-2B\£\/D], for GUE, 
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FIG. 2. In (a) and (b) the statistical distribution of £02 = ^0^2 over the intervals I = (— 0.14tt, — 0.06-7r) and 
T — (0.02-7T, 0.107r) respectively are compared with the prediction of the Gaussian hypothesis. The two intervals 
chosen correspond respectively to regions of phase angle which are scarred and (relatively) antiscarred. Figures (c) 
and (d) show the statistics of £22 = x\ in the same scarred and antiscarred intervals. Note the different scales in the 
scarred and antiscarred cases. 



When I = k we find in the GOE case that 

p(€kk' t 9) = 



i 



-1/2 



A- A 



y/2nC kk (0) 



-£kk/2Ckh(&) 



(16) 



which is just the Porter-Thomas distribution scaled so that (£kk) — Cfcfc(^), consistent with [Ol. 

In order to test these predictions with reasonable statistics, we create ensembles by selecting states with 
phase angles restricted to subsets X of the unit circle and compare their statistics with the accumulated 
distribution 



Pi{£,ik) 



1 



H(T) 



p(bk;O)d0, 



(17) 



where (i(I) = fjdd- These ensembles provide averages that are local in 8 while allowing sufficient data to 
make meaningful numerical tests. The resulting statistical distributions are compared in Fig. E3 with the 
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statistics of £02 = 2^2 and £22 = x\. In each case the statistics in two subintervals are shown: one covering 
a region where there is scarring (larger than average overlaps) and one in a region where there is antiscarring 
(smaller than average overlaps). In all cases the Gaussian hypothesis works well in describing the whole 
distribution. Note that similar agreement has been found for models with GUE statistics (which are not 
shown here in the interests of brevity) . 







FIG. 3. The distributions of xt, accumulated over the whole unit circle of phase angle, are shown in (a)-(d) for 
k — 0, 2, 4, 6 respectively. The insert in (a) shows the tip of the distribution. The Gaussian hypothesis (solid curve) is 
in good agreement with the quantum data (histogram) over the whole distribution. In each figure the dashed curve is 
the normal, unscarred, RMT prediction. The scarred distributions approach this RMT case as k increases (although 
not necessarily monotonically in general). 

We also show, in Fig. [j, distributions for the individual overlaps x^. The local distributions p(xk',0) 
are Gaussians scaled to that the variances are Ckk(&), analogously to (&). We then define accumulated 
distributions by integrating over an interval X with respect to 9, as in (|17|). In practice, the most natural 
ensemble is formed by averaging over the whole unit circle, corresponding to taking X = [0, 2tt). A sequence of 
such distributions is shown in Fig. |3| corresponding to fc = 0, 2, 4, 6. One sees that the distribution approaches 
the normal GOE prediction as fc increases but that for smaller values of fc the deviation from GOE is quite 
strong. Note that the Gaussian hypothesis describes quite accurately the statistics for moderate values of 
Xk as well as in the tails which correspond to the strong scarring limit. There is some underestimation in 
the antiscarring limit where Xk is small, but the overall agreement is good over the whole range of values of 
Xk- Similar or better agreement has been found in all the models we have examined. 

All the numerical evidence we have collected indicates strongly that the Gaussian hypothesis accurately 
describes the complete statistics of chaotic eigenstates in a neighbourhood of a periodic orbit. While mo- 
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tivation for extending the linear theory of scarring to a complete basis has initially come from a desire to 
treat tunnelling, we note that the results in this section are interesting from the point of view of scarring 
alone. The covariance matrix C(Q) provides a means of characterising the eigenstate statistics in the neigh- 
bourhood of a periodic orbit in a way that is independent of the choice of test state. We can accomodate 
the characterisation of statistics using a basis other than {\k)} simply by altering the covariance matrix 
using a similarity transformation. Ideally, any basis we use should in the interests of compactness cover the 
immediate neighbourhood of the periodic orbit (as measured by Wigner functions, for example) with as few 
members as possible, but we are otherwise then free to manipulate those basis elements at will. In addition 
to being essentially basis-independent in this way, such a characterisation is complete, in the sense if we 
choose a subset of the basis which covers a semiclassical neighbourhood of the periodic orbit, the statistics 
of any other measure of the eigenstate near the periodic orbit could in principle then be determined from 
the corresponding sub-block of the covariance matrix. 

Note finally that the deviations from RMT of the statistics of low-lying components of the eigenstate 
persist, and in fact remain fixed, in the semiclassical limit. This can be understood intuitively to derive from 
the fact that as Ti — > (or N — ► oo) the support of a harmonic oscillator basis state shrinks at the same rate 
as the semiclassical footprint of the periodic orbit. In E, the area occupied by each decays as 0(K) in the 
semiclassical limit. This semiclassical shrinking of the test states we use allows the scarring effect to remain 
fixed while other measures of scarring (such as the value of a wavefunction at a fixed point in space) lead to 
a vanishing effect in the semiclassical limit. 

IV. TUNNELLING STATISTICS WITH SCARRING 

The scarred wavefunction statistics described in the previous section are now used to derive distributions 
for the anomalous tunnelling rate statistics of [0,|3. The first step is to outline how the discussion for 
maps in the preceding section may be adapted to the statistics of states in potential- well problems and then 
we can describe how these ideas are used to construct distributions for tunnelling rates. 

While the formalism of tunnelling is presented in [n],|lJ| in terms of quantum maps and their eigenstates, 
quantisation of a potential-well problem leads more naturally to wavefunctions in a full Hilbert space and 
at first sight it is not obvious how to relate the two descriptions. Faced with a numerically-produced 
wavefunction, for example, there is no exact, purely quantum-mechanical way of defining a corresponding 
state in a finite-dimensional Hilbert space; all existing descriptions are semiclassical in nature. Our task is 
made easier by the fact that tunnelling depends only on properties of the state in a small region of phase 
space (or of a Poincare section in the reduced description) surrounding a tunnelling orbit. Within that region, 
the dynamical properties of the tunnelling orbit have been used in [|I5[ to give an explicit prescription for the 
definition of a Wigner function on a Poincare section. We will not repeat the details of that prescription here 
but state simply that a recipe exists for turning an eigenfunction ip n (x,y) of a Hamiltonian with potential 
V{x,y) into a Wigner function W n (y,p y ) on a Poincare section E defined by x = xq and on which (y,p y ) 
are canonical coordinates. This Wigner function is invariant near the tunnelling orbit under iteration of a 
semiclassical quantisation of the Poincare map [ |16| . With appropriate normalisation this construction allows 
us to compute overlaps, using 

i|fc) 2 = f dydp y W n (y, Py )W k (y, Py ), (18) 

between a reduced chaotic state \ri) and an eigenstate |fc) of a harmonic oscillator Hamiltonian h(y,p y ) 
defined on the Poincare section E. 

We now check that these overlaps are described statistically by the distributions derived in the preceding 
section. Until now the discussion has been given entirely in terms of simple quantum maps for which the 
chaotic spectrum is described in terms of eigenangles on the unit circle. In physical problems the spectrum 
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is more naturally described in terms of energy levels on the real line and we must connect the two pictures 
before proceeding. A detailed connection can be made using the transfer-operator approach of |16|, but for 
present purposes it is sufficient to identify the part of the spectrum corresponding to a narrow segment of the 
spectral envelope (so that we may partition states into ensembles with simple Gaussian overlap statistics) . A 
continuous-time analog of the phase angle 9 is provided by the action S(E) of the periodic orbit responsible 
for scarring fll3|] . We therefore partition the quantum states into ensembles for which 9 — S{E)/Ti lies in 
narrow intervals, mod 2-k. We assert that the statistics of the overlaps in (O) for each such ensemble should 
be a nonisotropic Gaussian such as described in (ffT 



P(Xq> 0.5 - V 




p(y) 




p(y) 




so 
o 



1.5 



1.0 - 



(d) 




FIG. 4. In (a) the theoretical distribution (solid curve) for overlaps with a harmonic-oscillator ground state is 
compared with quantum data (histogram) for the (y-symmetric) potential with (fi,u,a) — (0.1,0,0). The energy 
is fixed at the value E — 0.9 (where the potential saddle has energy E = 1) and ft is quantised in the quantum 
calculation. The histogram represents 392 states with 48 < 1/h < 75.77, consisting of 12 periods of the spectral 
envelop function. In (b) we show the corresponding distribution of rescaled tunnelling rates. In addition to the 
standard RMT prediction without scarring (short-dashed curve) we also show the Porter- Thomas distribution as a 
long-dashed curve. In (c) and (d) we present tunnelling rate statistics, using linear and log scales respectively, for 
the potential with (fj,, u, a) = (0.25, 0.40, 0.254). In this case states 1480 states with 40 < 1/h < 106.47, representing 
54 periods of the spectral envelop function, are used. This potential is not symmetric in y but the real extension of 
the tunnelling orbit is still periodic (inset). As in (a) and (b), the short-dashed curve represents a RMT calculation 
without scarring, the solid curve the scarring distribution and the long-dashed curve shows the Porter-Thomas 
distribution. 
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We test this assertion using the two-dimensional double-well potential 

V(x, y) = (x 2 — 1) + x 2 y 2 + /it/ 2 + vy + ax 2 y 

whose dynamics are predominantly chaotic near the saddle energy for appropriate parameter values and which 
we will later use to test tunnelling rate statistics. In practice we form statistical ensembles by fixing the energy 
and quantising h and in this way classical-dynamical parameters are held fixed throughout the spectrum. 
Our system does not allow the accumulation of sufficient data to test the distributions meaningfully for 
a fixed value of 9 = S(E)/Ti and we must instead pass straight to a distribution in which all states are 
combined, of the kind illustrated in Fig. 0. We show in Fig. 0(a) the accumulated distribution p(xq) of 
overlaps with a harmonic oscillator ground state (with k = 0) for the parameter values (fi, v, a) = (0.1, 0, 0). 
The potential in this case has a symmetry of reflection in y and therefore has a real periodic orbit connecting 
to a tunnelling orbit in the a;-axis. The harmonic oscillator used is the one which generates the monodromy 
matrix W of the tunnelling orbit and from which the tunnelling operator is constructed. The distribution 
derived from the Gaussian hypothesis describes the quantum data well and indicates that it is a reasonable 
basis from which to tackle tunnclling-ratc statistics. 

To describe tunnelling-rate statistics, it is convenient to begin once again with an ensemble taken from a 
fixed part of the spectral envelope (with 9 — S(E)/h fixed mod 2tt). Repeating the derivation of (jig), with 
Bo- relaced by the matrix T representing the scaled tunnelling operator, we obtain 

P(y,0) = ^- I™ e- l ™p{q,6)&q, 



where 



2n 



det ( I - ^-TC(9) 



-0/2 



(19) 



Our numerical calculations are made for GOE systems, in which case (3=1. (Note that while the mechanics 
of the calculation are similar to those leading to (flq), the quantities under investigation are very different 
in each case.) This distribution describes an ensemble of tunnelling rates taken with a fixed value of 9. As 
was the case for overlap statistics, it is not possible with our present model to produce enough data to test 
the distribution for narrow ranges of 9 and instead we make a comparison in which states from all parts 
of the spectral envelope are combined. Theoretically, we produce a combined distribution describing such 
collective statistics by averaging p(y, 9) over < 9 < 2n. 

Comparison of the resulting distribution with quantum data for the potential V{x, y), with the same pa- 
rameter values as in Fig. 0(a), is shown in Fig. 0(6) as a solid curve. Agreement with the scarred distribution 
is good and there is marked deviation from the prediction of standard RMT (short-dashed curve). Also 
shown for comparison is the Porter-Thomas distribution [13] (long-dashed curve). 

Systems with symmetry such as that used for Figs. 0(a) and H(o) are expected generically to show scarred 
statistics of the kind described in this paper. Nonsymmetric systems may also have scarred statistics, 
however, if parameters are tuned so that the real extension of the tunnelling orbit is periodic. We show 
tunnelling rate statistics in Figs. Q(c) and ||(d) for one such potential, corresponding to the parameter values 



(fi,i>,a) = (0.25,0.40,0.254). Quantum data for this potential were used in 11 to illustrate that dynamics 
could affect tunnelling statistics and the present theory gives a quantitative explanation of the qualitative 
observations made there. In Fig. [|(c) we compare the scarred p(y) (solid curve) with the standard RMT 
prediction (short-dashed curve) and the Porter-Thomas distribution (long-dashed curve). 

Note that while the scarred distribution appears not to differ strongly from the Porter-Thomas distribution 
in Figs. [|(6) and [|(c) , the difference is significant when plotted on a logarithmic scale. This is done for the 
nonsymmetric potential in Fig. 0(d) (and omitted for the symmetric potential in the interests of space). We 
also note that in the case of the symmetric potential the first hundred or so states (omitted from Figs. |](a) 
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and §(6)) give a distribution that seems better described by Porter-Thomas, though the scarred distribution 
unambiguously describes better the more highly excited states corresponding. This semiclassical limit seems 
to be approached more rapidly for the nonsymetric potential. 

We note that while the detailed discussion has been presented here for two-dimensional systems, the 



general theory will also work in higher dimensions. In particular, for example, (19) holds except that the 
matrices T and C(9) have additional structure. Generalisation of T to higher dimensions is straighforward 
(see pB[) and, for C(9), the generation of the autocorrelation function as described in the appendices may 
also be extended, though there are now more than the three classical parameters (TrVF, TrMo, Q), controlling 
quantum recurrence. 

The success of scarring theory in describing the tunnelling-rate distribution therefore confirms the hypoth- 
esis of ]TH that dynamical details of such low-dimensional systems can dominate the statistics of tunnelling 
and be responsible for strong deviation from RMT. The deviations we have seen are expected to be common 
in systems with discrete symmetry, such as the double well potential we have treated in Figs, ^(a) and 
\Q(b). Deviation will also be seen in more generic systems, however, as parameters are tuned to make the 
real extension of the tunnelling orbit periodic flllfl . Such sensitive dependence of experimentally-accessible 
data on system parameters and dynamics might be useful, for example, as a means of probing the internal 
dynamics of relatively complex systems. 

V. CONCLUSION 

We have shown that wavefunction statistics in the neighbourhood of a fixed point or short periodic orbit 
are well described by a normal distribution which allows for correlations between wavefunction components 
along different basis elements. We have also shown that this wavefunction distribution can be used to 
describe scarring anomalies previously detected in tunnelling-rate statistics. 

Scarring effects lead to strong deviations from the standard distributions of RMT in the neighbourhood of 
short periodic orbits. To characterise these deviations fully, we examined the components of the eigenstates 
of a chaotic map in the eigenbasis of a harmonic oscillator centred on a periodic orbit. The low- lying states 
of such an eigenbasis are localised around the corresponding periodic orbit and overlaps with a chaotic 
wavefunction are strongly affected by scarring. In appropriate subintervals of the chaotic spectrum, we 
have proposed that the joint probability distribution describing the corresponding components of chaotic 
eigenstates remain normal as in standard RMT but has a nontrivial covariance matrix (whose elements 
are completely determined by linearised dynamics around the periodic orbit). This hypothesis suffices to 
describe accurately all the measures of deviation from RMT that we have examined. The collective statistics 
of chaotic wavefunctions taken from a complete chaotic spectrum are described by superimposing these 
distributions. Deviations from RMT remain strong (and are effectively reproduced by the theory) in these 
combined ensembles, which are more natural to work with from a practical point of view. 

We believe that examining the wavefunction statistics in a complete basis such as this provides a promising 
technique for understanding scarring in general terms. The description is essentially basis-independent and 
also completely determines the nature of wavefunction statistics near a periodic orbit. 

The assumption of normal deviation from RMT also allows the effect of scarring to be incorporated very 
naturally into the calculation of tunnelling rate statistics. The theory was tested on a two-dimensional 
double-well potential with chaotic dynamics and found to describe the statistics of its eigenfeuntions and 
energy-level splittings well. Scarring can account for quite strong deviation from standard RMT in this 
system. In this way we might use tunnelling and the statistics of tunnelling to probe sensitively the internal 
dynamics of unstable systems. 
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APPENDIX A: GENERATING THE CORRELATION FUNCTION 



We calculate the correlation function Aikit) defined in section III using a generating function 

G(w,z,t) = Tie wa f (z)Uf ln , (Al) 

where a is an annihilation operator for the harmonic oscillator h and 

f(z) = e {lnz) ' h/h , 

reduces to the tunnelling operator when z = e~ a °. Each of the three operators in this sequence is the 
quantisation of either a linear symplectic transformation or a phase-space translation, possibly complex. The 
unitary operator Uf in quantises a real symplectic matrix Mq which linearises the real unstable motion near 
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(o- The operator T(z), which is nonunitary, quantises the complex symplectic matrix W(z) defined in (|l4|). 
Finally, the operator e wa can be interpreted as a translation along the complex phase-space displacement 



«~Vi(! 



The combined operator therefore quantises the affine transformation 

C h-> W{z)MlC + S(. 



(A2) 



Because the transformation is affine, the trace formula may be used to evaluate the right hand side (Al) 
exactly. It is shown in appendix M that this gives 



G(w, z, i) = 



1 



: exp 



(-1)"*' 



2 A(z,t) 



t&nip(t) 



where 



A(z,t) =TiW{z)M*-2= (-l)^f m (t)z + ^-^] 



(A3) 



and m(t) = e 1 ^ 1 ^ secip{t) as described in Sec. III. Note that G(w, z, t) is independent of %. Note also that 
since A(z, t) is complex, care must be taken in deciding the branch of the square root in ( |A3| ) — the branch 
is chosen using the Maslov index (computed from real dynamics), so that A(z,t) reduces to the standard 
formula for real maps when z — > 1. 



Expanding the right hand side of (A3) gives 

G(w,z,t) = ^(-l)"' it ^ T ^z fe+2 " +1 / 2 sin>yooIv;C fe " +1/2 (cosV) e <("+*+i/2)(*-M**) ) 



n=0 



fc=0 



where C^ (cos ip) are Gegenbauer polynomials defined by the generating function 



£cT 1/2 (*)* 



(1 - 2XZ + 2)n+l/2 ^ 

Comparing this with the expansion 

Co oo m 

G( W ,z,t) =££ ^_ zfe+ ™ + i/2 Ck\a m UL\~k) 

l — ' £ — ' rrj.l 



fc=0 m=0 



(A4) 



obtained by evaluating the trace (Al) in an eigenbasis of h leads us to the result 



(~k\a 2n Uf hl \~k) = ^jSin n t(;^^pC^ +1/2 (cos^) e *(k+V2)W-M*0+in*+*»r/2 



from which we determine AJ5: 2n fc(*) as gi ven in CP~3|) - Note that we find ^4}^ n (t) = unless I and fc are both 
even or both odd. This is because the linearised system enjoys an inversion symmetry even when the original 
system does not. 
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APPENDIX B: TRACING THE AFFINE TRANSFORMATION 



In this appendix we derive (A3) using the trace formula | 
mation 



We consider first the general affine transfer- 



(Bl) 



where M and Z = (a, b) are an arbitrary symplectic matrix and phase-space displacement respectively. For 
convenience we may initially suppose that M and Z are real and generalise to complex transformations later 
using analytic continuation. 



The transformation in (Bl) has a single fixed point 

C = -(m -iy x z. 

The trace formula then gives [E2| 

TiT(Z)U(M) = ■ 

where T(Z) and U(M) are the operators which quantise £ 1— > C + -^ an d C l— * ^C respectively and So is the 
action of the fixed point. By evaluating the generating function for M at the fixed point one can show that 

So = -Abq-aP) 



n(c,z) 
n(z,(M-iy l z), 



where £l(u, v) = u T Jv is the symplectic form and £ = (q,p). 



We specialise to the trace in (Al) by substituting the S( defined in (A2) for Z and W(z)AIq for At, giving 

G(w, z, t) 
where 



y /TT:W{z)M t 



E(z,t) = n c 



1 



W(z)M* 



<c 



and £e 



(B2) 



1 /9 

Notice that, because <5£ scales as % ' , there is no ^-dependence in the final result. 
We are free to work in a coordinate system in which K — I and then 



Consider first the case 



W(z) = e 



Ml = (-!)*" 



i In z J 



= z(I + U) + -(I-iJ). 



cosh pt e a sinh pt 
e^sinhp^ coshpt 



(B3) 



where a and p are constants. Substitution in (B2) gives then 

E(z,t) = 2(-l) M * cosher sinh pi 



A(z,t)' 
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where A(z, t) was defined in appendix A. Explicit calculation of A(z, t) for this Mq gives 

Q = — sinh a. 

Using the identity 

sec ip(t) = cosh pt + sinh crsinh pt — 1 + cosh crsinh pt 

we then get tan ?/>(£) = cosher sinh pi and hence (|A3|). 

Any hyperbolic matrix can be put in the form ( p33| ) following an orthogonal transformation which ro- 
tates coordinates by an angle 6$. The matrix W{z) is unchanged by such a transformation. However, the 
annihilation operator is multiplied by a phase 



and the exponent of the generating function by 

E(z,t) -> e 2i6 °Z(z,t). 

In systems with time-reversal symmetry we can demand that the eigenvectors |fc) be real and this allows 
us to fix the possible rotation angles 9q to multiples of ir/2. Among these, only 8q — and do — ir/2 give 
distinct results for E(z,t). The formula we have given applies to the case 0q = 0. In the case 0q = 7r/2, the 
formula holds as given if we rephase the eigenvectors according to \k) — » (— l) fc |A:). 

We can similarly use a convention for \k) such that ([1^) holds in the GUE case. This amounts in particular 
to adopting a convention in which the covariance matrix Cik{E) = (xix^) is real even though individual values 
of x\x\ are complex. This reflects the fact that even in systems without global time-reversal symmetry, the 
linearised dynamics from which A^£{t) is calculated is time-reversal invariant and we are free to choose the 
eigenstates |fc) in a way which reflects this. This is in direct analogy to the spatial symmetry mentioned 
in the previous appendix which the linearised dynamics also possesses even in the absence of an equivalent 
symmetry for the global dynamics. 
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